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Abstract 

The exact general solution to the Einstein equations in a homo- 
geneous Universe with a full causal viscous fluid source for the bulk 
viscosity index m = 1/2 is found. We have investigated the asymptotic 
stability of Friedmann and de Sitter solutions, the former is stable for 
m > 1/2 and the latter for m < 1/2. The comparison with results of 
the truncated theory is made. For m = 1/2, it was found that families 
of solutions with extrema no longer remain in the full case, and they 
are replaced by asymptotically Minkowski evolutions. These solutions 
are monotonic. 
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1 Introduction 



It is believed that quantum effects played a fundamental role in the early Uni- 
verse. For instance, vacuum polarisation and particle production arise from 
a quantum description of matter. However it is known that both of them can 
be modeled in terms of a classical bulk viscosity [Q. Other processes capa- 
ble of producing important dissipative stresses include interactions between 
matter and radiation |2| , quarks and gluons || , different components of dark 
matter ||, and those mediated by massive particles. It happens also due to 
the decay of massive superstrings modes ||, gravitational string production 
M [0 and phase transitions. 

Cosmological models with a viscous fluid have been studied by several au- 
thors. Some of the interesting subjects addressed by them were the effects of 
viscous stresses on the avoidance of the initial singularity |8j , the dissipation 
of a primordial anisotropy || , the production of entropy H , inflation and de- 
flation [|10| . The evolution of a homogeneous isotropic spatially- flat universe 
filled with a causal viscous fluid, whose bulk viscosity coefficient is related 
to the energy density by the power law ( ~ p m , and a truncated version of 
the transport equation for the viscous pressure, has been investigated thor- 
oughly |TT|| [FJI [131]. This model was based in the relativistic second-order 
theory of non-equilibrium thermodynamics developed by Israel and Stew- 
ard [Q. Their formulation employs fiducial local equilibrium magnitudes 
like thermodynamical pressure and temperature, related by an equilibrium 
Gibbs equation, but generalizes the expression of the entropy flux by the 
inclusion of quadratic terms in dissipative non-equilibrium magnitudes like 
bulk-viscosity pressure and heat flux. An alternative formulation of this the- 
ory, called Extended Irreversible Thermodynamics has been made by Pavon, 



Jou and co-workers |15j ||16|| . Here, non-equilibrium magnitudes like temper- 
ature and pressure are introduced instead, and a generalized Gibbs equation 
including the dissipative magnitudes is employed. The equivalence between 
both formalisms up to second-order departure from equilibrium has been 



demonstrated Jl7j 



Recently, it was shown that the truncated version of the transport equa- 
tion leads to pathological behaviour in the late universe, as the thermody- 
namical temperature increases with the expansion, leading to an unphysical 
heating of the Universe ]Tj|. Then, some authors have begun to employ the 
full version of the transport equation for the bulk viscous pressure, and made 
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comparisons with the results obtained with the truncated one [HJ [20J 
IS] [22] [24 1 [25 1 [26]. One important issue that arises in this context is 
the choice of the equations of state. For a dissipative Boltzmann gas it was 
shown that inflationary solutions appear not to arise, and it was conjectured 
that the inflationary solutions of the truncated causal theory are spurious 
and an artificial consequence of truncation [19|. In [^7| we have adopted 
also the full version of the transport equation but we have taken the simple 

and we have classified 



li 



forms of the equations of state used in |2jj and 
the possible behaviors near a singularity. In the present paper we complete 
the analysis of this cosmological model. We present the set of equations that 
describe the model in section 0. In section ^ we solve them when m = 1/2 
and give a detailed analysis of their solutions, including an analysis of the 
energy conditions |28[ . In section [| we study of stability of those asymptotic 
solutions that appear for m ^ 1/2 in the truncated causal and noncausal 
theories by means of the Lyapunov method. The conclusions are stated in 
section [|. 



2 The Model 

In the case of the homogeneous, isotropic, spatially flat Robertson- Walker 
metric 



ds 2 = dt 2 — a 2 {t)(dx\ + dx\ + dx\ 



(1) 



only the bulk viscosity needs to be considered. Thus we replace in the Ein- 
stein equations the equilibrium pressure p by an effective pressure @ 



H 2 



1 

3 P 



1 



H + 3H 2 = -(p-p-a 



(2) 



(3) 



where H— a/a, ' = d/dt, p is the energy density, a is the viscous pressure, 
and we use units c = 8nG = 1. As equation of state we take 



V = (l - i)p 



(4) 



3 



with a constant adiabatic index < 7 < 2, and for a we assume the transport 
equation, 

a + T& = -3(H - \ra (3H + I - | - ^ J (5) 

which arises from the causal irreversible thermodynamics theory, as the sim- 
plest way (linear in a) to guarantee positive entropy production |H| . 
Here £ is the bulk viscous coefficient, and following [[^J 

C = «p m (6) 

where a > and m are constants. Also, r is the bulk relaxation time and 
we choose the law 

r = i (7) 



which ensures causal propagation of viscous signals [j30fl . Finally, T is the 
equilibrium temperature and following [21j we choose 

T = np r (8) 

where k > and r > are constants, which is the simplest way to guarantee 
a positive heat capacity. In the context of standard irreversible thermody- 
namics, p, p, T and the number density n are equilibrium magnitudes, which 
are related by equations of state of the form p = p(T,n) and p = p(T,n). 
Further the thermodynamic relation holds [ffl 



' dp \ p + p T ( dp\ 
dn) T n n\dT) n ' 



(9) 



which follows from the requirement that the entropy is a state function. In 
our model, this relation imposes the constrain 

T - 1 

r = (10) 

7 

so that < r < 1/2 for 1 < 7 < 2. Using ©-©, we obtain 



7t2 o ol—m 

H - (1 + r)=- + - [1 + (1 - r) 7] HH + H 2 \H\~ 2m H- 
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q o2-m 

-( 7 -2)# 3 + - 1 H i \H\- 2m = U (11) 

for m^l/2, and 

HH - (l+r)H 2 + AH 2 H - BH 4 = (12) 

for m = 1/2, where A = 3/2 [1 + (1 - r) 7 + 2/70], B = 9/4 (2 - 7 - 27/70), 
and 70 = \^3a sgnH. It is easy to check that solutions of fli~2"D never change 
sign, and we will consider in the following only expanding evolutions, that is 
H > 0. Thus, taking into account (fLOf), we note that A > 0. 



3 The casern =1/2 

For the integration of (0), provided that r 7^ 0, it is convenient to make the 
change of variable H = y^ 1 ^. It turns this equation into 

y + Ay'ry + Bry 1 ^ = (13) 
which is a particular instance of the class of equations 

V + aif(y)y + 0C2f{y) J dyf{y) + a 3 f(y) = (14) 
that linearises through the transformation [13| |27| |3^] 

z = J dyf{y), rj = J dtf{y) (15) 

becoming 

cP z dz 

drf + ai dr] + + " 3 = ( 16 ) 

Thus, for «! = 1, «2 = (r — 1)B/A 2 and 03 = we obtain the general solution 
of flT2|) in parametric form 

#(77) = [i?z(r])]^ 
A*fa) = j J dvHirj)- 1 (17) 
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where At = t — to, and to is an arbitrary constant and R is a constant. For 
the case r = 0, that is isothermal matter, it is easy to verify that flT7|) is also 
the general solution. The constrain ( |XD[ ) implies that 

z(r]) = C + exp (X + rj) + C_ exp (A_r/) (18) 

1 /2 

where X± = (1/2) [— 1±(1 — 4«2) ] are real, and C± are arbitrary integration 
constants. When any one of them vanishes, we obtain the one-parameter 
families of solutions, which have explicit form: 



#4 



H±(t) - 
- vo/At, 



u ± /At, 



U±= 2B 



-A± 



A — 4B (r — 1 



B^O 
Ho, B 

1/2" 



;y 



1 — r 
A 



(19) 
(20) 

(21) 



and Ho is an arbitrary positive constant. As < v + < (2/3) (v^ — 1) and 
l/6<t / o<2/9, the behavior of a + (t) is Friedmann. On the other hand, the 
de Sitter behavior for B = has already been considered in |18[ |22 . 



3.1 Two-parameter families of solutions 

As shown in appendix A, for some values of 7 and 70, a(t) can be written 
in closed form in terms of known functions. However, for arbitrary values 
of 7 and 70, we need to study the solution in the parametric form flTT|). 
Depending on the case, the behavior of H(t) for small or large time is given 
by the behavior of z(r]) near a zero point or for ±00. In the first case, z ~ rj 
for r] -»• 0, H ~ l/(At) 1 / r , provided r ^ 0, or H ~ exp(HAt) if r = 0. Thus 



a(t) ~ aoexp [KAt r 



r > 



a(t) ~ ao exp 



-exp (KAt) 







(22) 



where ao and are arbitrary integration constants. In the second case, either 
the first or the second term in (|TB| ) dominates, and the leading behavior of 
H(t) is given by (|19|) (or (|20|), if B = 0). Thus, we obtain the following 
classification of the two-parameter families of solutions: 
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A. The evolution begins at a singularity with a Friedmann leading be- 
havior as ([i9|+) for B ^ or (f20|+) for B = 0, and so there are particle 
horizons. Then the scale factor expands with an asymptotically Minkowski 
behavior for large time like (|22|). 

B. The evolution is asymptotically Minkowski in the far past with be- 
havior (^), then it expands and its behavior in the future is either: 

1. asymptotically Friedmann, as (|19| -) for B < 0; 

2. asymptotically de Sitter, as (p0|-) for B = 0; 

3. divergent at finite time, with leading behavior (|19|-) for B > 0. 

C. The evolution begins at a singularity with a Friedmann leading be- 
havior as ([i9|+) for B ^ or (|20|-|-) for B = 0; and so there are particle 
horizons. Then the scale factor expands, and its behavior in the future is like 
B. 



3.2 Energy conditions 



Let us consider whether this viscous fluid obeys the energy conditions when 
m = 1/2 f28|. Using (H)©, the dominant energy condition (DEC) p > 



\p + er| implies —3H 2 < H < 0, while the strong energy condition (SEC) 
p + 3p + 3a > implies H + H 2 < 0. We find that DEC is violated part 
of the time in the following families: Bl, in the far past as well as for large 
times (if V- < 1/3); near the singularity in CI, C2 and C3; C3, near the 
"explosion". DEC is violated always in families A, B2 and B3. SEC is 
satisfied always in families A and Cl (if v_ < 1). It is violated part of the 
time in the families: Bl, in the far past (if z/_ < 1); for large times in Cl 
(if V- > 1), C2 and C3. SEC is violated always in families Bl (if z/_ > 1), 
B2 and B3. Thus, we find that all these two-parameter solutions violate the 
energy conditions sometime. 



4 The case m^l/2 



We will make use of the method of the Lyapunov function [34] to investigate 
the asymptotically stability of the de Sitter and asymptotically Friedmann 
solutions that occur in the noncausal and truncated causal models. 
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4.1 Stability of the de Sitter solution 

Equation (|TT| ) admits a de Sitter solution for 7 < 2 [[IS 



Ha 



o a 



' 2-7 * 
, 2 7 . 



(23) 



To study its stability, if r > 0, we make first the change of variable H = y 1 ' r 
in (0) 



3 1 gl— m 

£ + - [1 + (1 - r) 7] y'-ry + ^ (m ~ 1} ?/+ 

2 a 



9 

- (2 - 7) n/ 



1-2 3 2 ' m 7 1. 



3 1 2m 

7 + — = 



and then we rewrite it as 



(24) 



d_ 
~dt 



\v 2 + V(y) 



yl—m 



-[l + (l-r) 7 ]y " + 



2(m-l) 



■2 



(25) 



where 



^(2/) = "7777 TV"?/ 



32-m 2 



2r-3+2m 



8(r- 1 
for 2r — 3 + 2m 7^ 0, and 



2a(2r-3 + 2m)' 



rr , . 9(2 — 7) r 2 2(r-i) 32 +r 7r 1 y 
8 (r - 1) 2a y 



(26) 



(27) 



for 2r — 3 + 2m = 0. We see that this potential has a unique minimum for 
y > at 2/0 = -^o~ r provided m < 1/2, while it has a maximum for m > 1/2. 
Also, taking into account (pf) , the right hand side of (|25|) is negative definite. 

On the other hand, if r = 0, we make the change of variable dt = dx/H 
in (0), and we obtain 



d_ 

dx 



~#' 2 + V(H) 



-H 



ol— Tn 

3 + E l - 2m 

a 



(28) 



where ' = d/dx, and 
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}2— m 



V(H) 



-H 2 + 



2a (3 - 2m) 



H 



3-2m 



m ^ 



^(F) = --iy 2 + ^ln— , 
V ; 8 2a # 



3 

m = — 
2 



(29) 



(30) 



This potential has a unique minimum for H > at ifo provided m < 1/2, 
while it has a maximum for m > 1/2. Also the right hand side of ( PE| ) is 
negative definite. 

Thus we find that an exponential inflationary regime is asymptotically 
stable for t — > oo for any initial condition H > provided that m < 1/2, 
but this regime becomes unstable for m > 1/2. This result improves over 



previous studies based in small perturbations about solution ([23]) |I8| |22 



4.2 Stability of the asymptotically Friedmann solution 



For m > 1/2 it is easy to check that ( |iT|) admits a solution whose leading 
term is 2/{?> r )t). To study its stability, if r > 0, we make the change of 
variables 



H 



[«(*)]' 



2m- 1 



in (|Tl|) which takes the form 



d 

dz 



l -u ,2 + U(u )Z ) 



D(u, u' ', z) 



(31) 



(32) 



where ' = d/dz. We consider that u lies in a neighbourhood of u = 
(37/2) _1 ^ r . Thus, when z — > 00 



gl— m^2 



2a(2m- l) 2 
for 2r + 2m ^ 2, 3, 



2(r+m-l) 



37M 



2r + 2m-3 



r + m — 1 2r + 2m — 3 



1 ^f 1 

- + °(- 

z Vz 



(33) 
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and 



D (u, u' z) 



;1— m 



2(m-l) 



a(2m 



-u 



u' 2 + O - 



On the other hand, if r = 0, we make the change of variables 

exp [u(z)} 



H 



z = t 



2m- 1 



(36) 



(37) 



in (11) which takes again the form (^), but now we consider that u lies in a 
neighbourhood of Uq = ln(2/3), and when z — ► oo 



gl— 771 

U(u,z) = o 

V ' 2a(2m-l) 2 



2 e (3-2m)« e (2-2m)u 



3 - 2m 



m 



2 Z A 



for 2m ^ 2, 3, 



U(u,z) 
U(u,z) = 



1 /3 



a V2 



m 



4^ 



-« + e u - + 0(— ) 



2m 



and 



D(u, u', z) 



> 1—771 



_ e (2-2,n) V 2 + Q J_ 



(38) 

(39) 
(40) 

(41) 



a(2m - 1) 

As U(u, z) has a unique minimum at uq for any m ^ 1/2, and D(u, u', z) is 
negative definite for m > 1/2, we find that solutions with leading Friedmann 
behaviour a ~ t 2 ^ 3 ^ when t — > oo are asymptotically stable for m > 1/2 , 
but become unstable for m < 1/2. 
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5 Conclusions 



We have improved our previous analysis of a cosmological model with a 
causal viscous fluid by considering the full transport equation for the bulk 



viscosity instead of the truncated one fl3| . We have kept the equation of 
state that relates the equilibrium pressure and energy density, and we have 
added a power-law relationship between the equilibrium temperature and 
energy density. 

When m = 1/2, the splitting of the large time asymptotic behavior of 
solutions in terms of sgn(— B) follows closely that of the truncated model 
in terms of sgn (7 — 70 ), which in its turn resembles the classification for 
the noncausal solutions. However, families of solutions with extrema no 
longer remain in the full case, and we find instead asymptotically Minkowski 
evolutions. Moreover, all solutions are monotonic. 

The singular behavior for m = 1/2 is also very similar to that of the 
truncated model because, excluding fll^-), singular solutions have particle 
horizons. The only difference is that u + in fl2l|+) is smaller than the cor- 
responding magnitude in |T3. Similar result has already been verified for 



m > 1/2 

As in the truncated model, we demonstrate that there is no value of 
m for which there is a stable expanding de Sitter period in the far past. 
This further support the conclusion of that causality avoids the deflationary 
behavior proposed by Barrow [[?]. 

If m < 1/2, a stable exponential inflationary phase occurs in the far 
future for 1 < 7 < 2; the condition B = is required if m = 1/2, and such 
a behavior is unstable for m > 1/2. This is the same structural behavior 
already found in the causal truncated and noncausal models. It has been 
noted however that the inflation rate in the full causal theory is different, 
lower for m < 1/2 and greater for m > 1/2. Also, there is no stiff-fluid 
viscous exponential inflation, unlike in the truncated or non-causal theories 

It has been conjectured that viscosity-driven exponential inflationary so- 
lutions are spurious and an artificial consequence of use of causal truncated or 
noncausal theories ]nj. This suggestion arises from study of a cosmological 
model containing a dissipative Boltzmann gas. Our results show however that 
this kind of conclusions does not generalize to other models where the vis- 
cous fluid is described by a different equation of state. Moreover, they rather 
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suggest that some features of the solutions, like their asymptotic behavior 
depend more strongly on the equation of state than on the thermodynamical 
theory employed. 

If m > 1/2, we find that the perfect fluid behavior a ~ t 2 /^) for t — > oo is 
asymptotically stable. This arises because the viscous pressure decays faster 
than the thermodynamical pressure. However, if m — 1/2 and B < 0, both 
pressures decay asymptotically as t~ 2 and the exponent becomes v_. The 
perfect fluid behavior becomes unstable if m < 1/2. Here also, no structural 
change occurs with respect to the truncated causal or noncausal models. 



6 Appendix A 



The equation ( |i"2"D has a two-dimensional Lie group of point symmetries, 
whose generators are d/dt and td/dt — Hd/dH. We look for a simple trans- 
formation to a new dependent variable v that is linear in v, preserves this 
group of symmetries and keeps the order of the differential equation. This 
technique can be extended to some polynomial autonomous differential equa- 
tions arising when we consider quantum effects due to vacuum polarisation 
terms |35|. Thus, we are led to the transformation 



H = fH 2 + gH- 



(42) 



where / and g are two constants determined by the requirement that the 
differential equation in v(t) is linear. Inserting fl42[) in (]TJ) we obtain 



l-r) f + Af-B\ H 4 -g[(2r-l)f-A] -H 3 - 
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— + ( gr + r 



v 2 



H 2 = 



(43) 



Then, provided we choose / = A/ (2r — 1), g — — 1/r, and we impose the 
constrain 



Br 
'A 2 



(44) 
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equation (O) turns into v = 0. To integrate (f42|), we make the change of 
variable H = —iij (fu), an it becomes 



u lv 
- + — 

u r v 



which has the first integral uv 1 ' r 
constant C. Thus we obtain 



= (45) 
C, with some arbitrary integration 



and the scale factor 



H 



a(t) = do 



2r) 



_ i 

V r 



A J dtv~ 



I Ail 



K 



(46) 



(47) 



where v — (1 — 2r)/A and K is an arbitrary integration constant. We 
note that, using (|i~0|), ( fi4"D becomes a curve on the parameter plane (7,70). 
This explicit general solution can also be derived from a more comprehensive 
mathematical framework, as shown in . 

For B = 0, an implicit expression for the solution of (|12D exists. In effect, 
equation (|T^) admits a first integral 



Ar r=i 

y-\ y r 

r — 1 



E 



and we obtain 



1 /iV r 



At = — [ — 
H \H 



1F1 1, 



2r - 1 Ar 



H 



l- r \ 



At 



r / a \ i-r 



AH \a Q 



iF 1 



r — 1 ' r — 1 ' r — 1 \Hq 
I 1 2-r /a 



1 — r'l — r'l — r'Vao 



(48) 

(49) 
(50) 



where Hq, ao are arbitrary integration constants. 
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Figure Captions 



Fig. 1. Selected solutions for 7 = 1.1 (B > 0). The curves of class A 
and C3 correspond to (f47[) with K = —1. The curve of class B3 corresponds 
to (0) with K = 1. 
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